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Abstract 

I review progress made in understanding the nature of the broad-line region (BLR) of active galactic nuclei (AGNs) 
and the role BLRs play in the AGN phenomenon. The high equivalent widths of the lines imply a high BLR covering 
factor, and the absence of clear evidence for absorption by the BLR means that the BLR has a flattened distribution 
and that we always view it near pole-on. The BLR gas is strongly self- shielding near the equatorial plane. Velocity- 
resolved reverberation mapping has long strongly excluded significant outflow of the BLR and shows instead that the 
predominant motions are Keplerian with large turbulence and a significant net inflow. The rotation and turbulence 
are consistent with the inferred geometry. The blueshifting of high-ionization lines is a consequence of scattering off 
inflowing material rather than the result of an outflowing wind. The rate of inflow of the BLR is sufficient to provide 
the accretion rate needed to power the AGN. Because the motions of the BLR are gravitationally dominated, and the 
BLR structure is very similar in most AGNs, consistent black hole masses can be determined. The good correlation 
between these estimates and masses predicted from the bulge luminosities of host galaxies provides strong support 
for the similarity of AGN continuum shapes and the correctness of the BLR picture presented. It is concluded that 
although many mysteries remain about the details of how AGNs work, a general overall picture of the torus and BLR 
is becoming clear. 
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1. Introduction 

A broad-line region (BLR) is present in all AGNs ac- 
creting at moderate- to high-Eddington ratios. BLRs are 
important both because they are our best probe of how 
AGNs work and because of their potential for readily 
providing masses of supermassive black holes (SMBHs) 
back to the earliest times of galaxy formation. How- 
ever, in order to be able to use BLRs to reliably es- 
timate the masses of SMBHs it is essential to under- 
stand the structure and kinematics of BLRs. Over the 
last four decades there have been wide-ranging and, not 
infrequently, mutually contr adictory views of the nature 
of the BLR (see reviews bv | Mathews & Capriottilll 985. 
Osterbrock & Mathews 1986, and Sulentic et al. 2000). 
However, I believe that the situation is improving. I re- 
view here what I consider to be the clearest pointers to 
the underlying structure and kinematics of the BLR and 
I argue that, while there are certainly many interesting 
problems remaining, the basic picture is now becoming 
fairly secure. I furthermore believe that this picture ap- 
plies to all BLRs because BLR equivalent widths and 
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line ratios are remarkably similar, especially in the ul- 
traviolet. 



2. The structure of the broad-line region and torus 

The two most basic questions about the BLR are 
"what does it look like?" and "how is it moving?" The 
traditional picture of the BLR of an AGN for over 40 
years (and one which is widely depicted in cartoons of 
AGNs) has been that there is a central source emitting 
ionizing radiation roughly spherically, and that it is sur- 
rounded by a roughly spherical mist of cloudlets. This 
is depicted in the left-hand-panel of Fig. 1. Each indi- 
vidual cloud, if it is big enough, will have a structure as 
shown in the right-hand-panel of Fig. 1 . It will be highly 
ionized on the front, and if it has a high-enough column 
density, it will be mostly neutral on the back. The front 
emits high-ionization lines such as Hell, He I, OVI, 
N V, and C IV, while the back emits low-ionization lines 
such as Mgll, Call, 01, and Fell. All these lines are 
well-known in AGNs. 
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Figure 1: The left frame shows a cartoon of a common traditional 
view of the BLR. The right frame shows a schematic close-up of an 
individual cloudlet. 
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Figure 2: Relative emissivity of four lines in a typical BLR cloudlet 
(normalized to a maximum emissivity of one for each line) versus 
distance (in arbitrary units) from the ionized face of the cloud.. 



The emissivity of each line as a function of distance 
from the front of the cloud can be calculated with the 



photoionization code CLOUDY (iFerland et all H998). 



Fi g. 2 shows emissivitie s for some well-known lines. 



Baldwin et al.l (1 19951) showed that the sum of contri- 
butions from clouds with a distribution of cloud proper- 
ties (densities and distances from the center) will auto- 
matically produce a total spectrum similar to what is ob- 
served from AGNs. This is the so- called LOC modelQ 
This was important because it showed that no "fine- 
tuning" of cloud conditions was needed to explain AGN 
spectra. 

Despite the success of the traditional picture in 
general, and the LOC model in particular, in explaining 
the overall spectrum of an AGN, the prob lem with this 
picture (see Gaskell. Klimek. & Nazarovall2007l) is that 
to explain the strengths of the BLR lines the covering 



factor has to be large (50% or so), yet if this is so, and if 
the cloudlets are covering the central source uniformly, 
we ought to see Lyman continuum absorption by the 
BLR clouds. In fact Lyman continuu m absorption due 



to the BLR is never convincingly seen ( Antonucci et al 



1989 - see discussion in MacAlpine 2003 and 
iGaskell. Klimek. & Nazaroval 12007). We believe, as 



proposed by S. Phinney (see lAntonucci et aT 1989) 



that the need for a high covering factor plus the lack of 
Lyman continuum absorption requires the BLR to have 
a flattened distribution and requires us to be viewing 
it through a hole. This conclusion is supported by re- 
covery of what is called the "transfer function" of some 
lines (the transfer function is the temporal response of a 
line to a delta- function event in continuum light curve). 
Transfer functions for low-ionization lines have always 
implie d that there i s little or no gas along the line of 
sight dKrolik et all Il99ll: iHorne. Welsh. & Peterson! . 

Mannucci. Salvati. & Strangal 1 19921: 
Pijpers & Wanders, Il994h . and thus that at least 
the low-ionization gas in the BLR has a flattened 
distribution. 

Having a high overall covering factor but a flattened 
distribution means that near the equatorial plane there 
will be a close to a 100% chance that any path will 
intersect a BLR cloud. The clouds will thus be self- 
shielding. Radiation from the central source can freely 
escape near the axis of symmetry, but is strongly diluted 
in the equatorial plane. This is schematically illustrated 
in Fig. 3. 

It is easy to calculate the average radial dependence 
of the ionization and the emissivities of all the lines 
coming from cloudlets with a distribution such as in 
Fig. 3. The ionization structure of a single cloud in 
CLOUDY is now spread out in radius as illustrated 
schematically in Fig. 4. The horizontal axis in Fig. 2 can 
now be read as distance into the BLR rather than dis- 
tance into an individual cloud. Our model is in fact very 
simila r to the old "filling factor" model of Mac Alpine 
dl972h . 

Th e earliest reverberation mapping of multiple 
lines dGaskell & Sparkei Il986h showed that the high- 
ionization lines were coming from smaller radii than 
the low-ionization lines. High - ionization lines were 



also wider (e.g., IShuderi 11982c Mathews & Wampler 
19851) . The radial ionization stratification of the 
BLR has been well confirmed by later reverberation 
mapping. The best reverberation-mapped AGN is 
NGC 5548. The horizontal axis of Fi g. 5 (taken 
from Gaskell. Goosmann. & Klimekl 120081) shows the 
reverberation-mapping time lags (i.e., the effective 



Ostensibly from "locally optimally emitting clouds." 



radii) for lines of a variety of ions from IClavel et al 
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Figure 3: Schematic cross section of the BLR and torus in a plane 
through the axis of symmetry. The torus is on the right. Ionizing 
radiation is attenuated in the equa torial plane, but can freely e scape 
near the poles. Figure from Gaskell, Klimek, & Nazarova ( 2007) 
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Figure 4: Cartoon of the relationship between a traditional cloudlet 
model (top two thirds of the diagram) and the self-shielding model 
(bottom third). The different shadings symbolize three regions pro- 
ducing lines of differing degrees of ionization (cf. Fig. 3). . 
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Figure 5: Observed lags (jobs) versus predicted lags (r pre d) for 
NGC 5548. The open triangles are the predictions of the LOC model, 
an d the solid dots are the pred i ctions of the self-shielding model 
of IGaskell. Klimek. & Nazaroval d200 7|). The lines are OLS-bisector 
fits (jlsobe et all 1 1 99oh . Figure from iGaskell. Goosmann. & Klimekl 
(2008). 



dl99ll) . |Peterson et alJdl99lh. andlBottorff et alJd2002h . 
(See table in lGaskelL Klimek. & Nazarovall2007l for de- 
tails). The observed lags cover an order of magnitude 
in radius. NGC 5548 is not unique in this regard: an 
iden tical range of radii has also been found for Mrk 110 
by iKollatschnyl (120031) . The vertical axis shows the 
lags predicted for the same lines by the L OC model 
dKorista & Goadl bOOd; liottorff et all l2002h . It can be 
seen that while there is a correlation, these predicted 
lags cover a much smaller range of radii. The reason for 
this can be appreciated in Fig. 4. In the LOC model (top 
part of the figure) every cloud has a highly-ionized front 
part; the clouds just differ in the degree of ionization. 

The self- shield ing model 

(IGaskell. Klimek. & Nazaroval 120071) solves the 
problem of why the ionization stratification is so strong. 
In the self- shielded model (bottom of Fig. 4) there is a 
clear spatial separation of the differing ionizations. It 
can be seen in Fig. 5 that the self- shielding model gives 
g ood agreement with the observed lags. 

iNetzer & Lao j ( 1993 ) made the important sugges- 
tion that the outer edge of the BLR coincided with 
the dust sublimation radius of the torus. Reverbera- 
tion mapping observations show that the low-ionization 
gas in the BLR i ndeed extends out to the dusty torus 
(Sug anuma et all 120061 IGaskell. Klimek. & Nazaroval 
120071) . The covering factor of the torus can be cal- 
culated statistically from the ratio of type-1 (face- 
on) to type-2 (edge-on) AGNs, and directly for in- 



Figure 6: A schematic view of the BLR and torus of an AGN 
in a plane through the axis. The figure is approximately to scale 
(except that the black hole is shown too large.) Figure from 
Gaske ll. Klimek. & Nazarov3 d2007h . 



Figure 7: Computer-generated renditions of the appearance of the 
BLR and clumpy torus of an AGN seen from a "Seyfert 1.9" viewing 
position (lower right), and a cutaway view from the same position 
(upper left). The cut makes a 45 degree angle to the projection of the 
line of sight onto the equatorial plane. Computer renditions by Daniel 
Gaskell. 



dividual objects from the strength of the thermal 



emission (e.g., LM aiolino et al. 



2007). 



We argue 



(IGaskell. Klimek. & NazarovaL [2007) that the covering 
factors of the BLR and torus have to be the same. This is 
because if the torus has a lower covering factor than the 
BLR we would see the BLR in absorption against the 
central continuum source in some objects near the type- 
2 viewing position. This is never seen. On the other 
hand, if the BLR has a lower covering factor, some of 
the dusty torus will see direct radiation from the central 
source. This cannot be the case for much of the torus 
because it would then be unable to exist as close in as is 
seen. 

The overall picture we get of the torus and BLR is 
indicated schematically in the cartoon in Fig. 6 and in 
the computer generated renditions shown in Fig. 7. The 
best description of the appearance is to say that the 
BLR and torus look like a bird's nest. This picture is 
identical to that favored for totally independent reasons 
in an unfortunately almost totally ov erlooked paper by 
Mannucci. Salvati. & Strangal dl992k . They inferred a 
"bird's nest" geometry from a combined analysis of line 
profiles and transfer functions in NGC 5548. The posi- 
tions of masers in NGC 1068 also provides support for 



tions or masers m J Nul iUba also provide 
a thick BLR-torus dGreenhill et all 1 19961) . 



3. THE KINEMATICS OF THE BLR AND TORUS 

3.1. Determining the Direction of Motion 

The kinematics of the BLR have been a long-standing 
problem. It has been known from the earliest days of 
AGN studies that the lines are very broad (for a review 



of the earliest literature see Seyfertl 1943 ). but Doppler 
shifts only tell us the motion of gas along the line of 
sight. To know whether the gas is inflowing, outflowing, 
moving in random virialized orbits, or in more planar 
Keplerian orbits in a disc we need to know the line-of- 
sight velocity as a function of position relative to the 
black hole. 

The discovery of narrow intrinsic absorption in NGC 
4151 dMavallL Il934l: lAnderson & KraftL Il969h andof 
broa d absorption lines (BALs) in PHL 5200 (ILynds . 
1967b proved that some gas was outflowing from AGNs. 
However, BALs commonly extend to velocities sev- 
eral times higher than those observed for the BLR 



in th e same objects (see, for example, iTurnshek et al 
Il988k . so it is not clear that there is necessarily any 
connection between BALs and BLRs. The case for 
an outflowing BLR w a s stren gthen ed thoug h when 



an o utflowin g bLK w as strengthened tnougn wnen 
Blu menthal & Mathews! d 1975k and iBaldwinl d 1975k 



showed that a radiatively-accelerated outflow could re- 
produce t he observed line profiles well in some objects. 
However, Capriotti. Foltz. & Byardl dl980k showed that 
other models could provide comparably good fits to 
broad-line profiles, and so demonstrated that fits to indi- 
vidual line profiles alone could not uniquely determine 
the kinematics. 

More progress was made by c omparing lines of dif- 
fering ionizations. Gaskelll d 1982k discovered that high- 
ionization broad lines were blueshifted with respect to 
low-ionization lines, and pointed out that this requires 
there to be radial motions plus some source of opac- 
ity. Thi s blueshifting has now been widely confirmed. 



Gaskelll dl982k suggested that the blueshifting could 
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Figure 8: Velocity-resolved reverberation mapping. Because of light- 
travel-time effects, the gas on the near side of the AGN is seen to re- 
spond to continuum changes first. For the the hypothetical outflowing 
BLR illustrated here, the blue wing of a line would vary first. 



be the result of a "disk- wind" model where the high- 
ionization line s arise in a wind ou t flowin g above the ac- 
cretion disc. Wil kes & Carswelll dl982h pointed out a 
problem with any purely radial motion: the profiles of 
C IV and Lyman a were observed to be very similar, 
yet, for optically-thick clouds, Lyman a is emitted very 
anisotropically. To satisfy this constraint the clouds ei- 
ther had to be optically thin, or not moving purely radi- 
ally. 

Obviously the question of the direction of motion 
could be settled if it could be determined which gas was 
on which side of the black holes. The best way of doing 
this is through velocity-resolved reverberation mapping 
dGaskellL[l988h . How this works is illustrated in Fig. 8. 
Surprisin gly, velo c ity-re so lved reverberation mapping 
results 



lgly, velo c ity-re so lved r everbera tion n 
(Gask eil 1 19881: iKoratkar & GaskellL 



Crenshaw & Blackwef 
1991a ; bx; [korista 



kwell, 
a et all 
1996; ! 



1990; 
19951: 



198 



Koratkar & Gaskeil 
Done & K rolikL 1 19961 



Ulrich & Hornel Il996l: ISergeev et all Tl999) strongly 
ruled out significant outflow of both high- and low- 
ionization lines (see example in Fig. 9). 

Ruling out significant outflow of the BLR was impor- 
tant not just because of what it said about how AGNs 
work, but also because it meant that the BLR motions 
were gravitationally dominated. The BLR could thus 
be used for determining the masses of the central black 
holes. This permitted the first reverberation mapping 
determin ations o f black hole masses and Ed d ington 



ratios dGaskellL 1988 



iKoratl 
Il99(j 



Koratkar & Gaskeil, 1989; 
Koratkar & GaskellL 



Cren shaw & Blac kwellL 

1991allbllch l 2 l Note that while gravity dominates BLR 
motions, the simple fact that radiation pressure was at 
one time considered to be dri ving BLR motions (e.g., 
Blument hal & Mathews! 1 1975k should warn us that ra 



diation pressure might not be negligible (IMarconi et al 
2008). 



While the velocity-resolved reverberation mapping 



2 The BLR w as first used to estimate masses of AGN black holes 
bv lDibail fl977h who estimated BLR sizes using photoionization con- 
siderations. At that stage, o f course, there was no clear e vidence that 
the BLR was virialized. See Bochkarev & Gaskeil ( 2009). 
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Figure 9: The cross-correlation function for the blue and red wings of 
the Mg II line in NGC 4151 as a function of time delay. The predicted 
peak in the correlation function for pure outflow (blue wing varies 
first) is shown by the arrow. It can be seen instead that the strongest 
correlation is for near zero delay (what is expected for virialized or 
Keplerian motion), but with the red wind leading by a small b ut sig- 
nifican t amount thus implying some net inflow. Figure from Gaskeil 
Il988h . 



results were good news for the new black-hole-mass- 
determination industry, they created a problem for 
the generally accepted "disk- wind" explanation of 
the blueshiftings of high-ionization lines. Disk- wind 
models are very theoretically appealing, and strong 
blues hiftings have been taken as signs of strong winds 
(e.g.. iLeighly & Moord 120041) . However, at the same 
time, people working on black hole mass determina- 
tions were firmly believing that they were using viri- 
alized lines ! This has almost caused AGN observers to 
suffer from multiple-personality disorder HQ 

It is very difficult to finesse a disk-wind model to 
fit the velocity-resolved reverberation mapping con- 
straint, all the more so since outflow w as first ex cluded 
for the / z/ff/?-ionizatio n_ C I V line (GaskdJ n 
Koratkar & Gaskeil, 



1985 



1988; 



Crenshaw & Blackwell, 



19901: lKoratkar& GaskellL ll991aHblch . We believe, 
however, that there is a simple solution to the problem: 
the opacity needed to cause the blueshifting is not pri- 



3 or at the very least to fear that, like the White Queen in Alice 
in Wonderland, they might have to believe in six impossible things 
before night lunch! 

4 The narrow-line region is creating a similar problem. People who 
study extended narrow-line emission associate it with jets and out- 
flows, while other people are using narrow-li ne velocity widths as a 
proxy for the stellar velocity dispersion (see Gaskeil 2009a). 
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Figure 10: Cartoon illustrating why scattered photons are blueshifte 
when scattered off a reflector which is approaching the source of phc 
tons. The person on the right sees her reflection (far left) in the mirroi 
If the mirror is approaching her, then the image is approaching he 
twice as quickly. 




Figure 11: Cross sections in a plane through the axis of symmetry of 
the two scattering region geometries modeled in Fig. 12. 



marily absorption but scattering (IGaskell & Goosmann 



2008). Electron scattering had in fact been considered 



in the late 1960s t o be a significant source o f line broad- 
ening in AGNs dKaneko & QhtaniL 1 19681: IWevmannl 
1970l: iMathisL Il970h , but the idea fell out o f favor 



with the success of the Blumenthal & Mathews! ( 1975 ) 
radiative acceleration model in fitting pr ofiles. It has , 
however, long been well kn own (e.g., Edmonds! 1950l: 
Auer & van Blerkoml 19721) that scattering off regions 
with a net radial motion produces line shifts. For 
an infalling scattering medium, photons gain energy. 
This is explained in Fig. 10. The process is similar 
to the well-known Fermi acceleration process. The 
effect of scattering off radially moving materia l in 
AGNs was consid e red bylKallman & Krolik (1986) and 



Ferrara & Pietrini ( 1993). 



As can be seen in Fig. 9, velocity-resolved rever- 
beration mapping not only excludes outflow, but it 
also shows that there is a slight inflow. Initially this 
result only had -90% confidence for any one line in 
one object, but it has been found repeatedly for many 
lin es in many objects now a nd thus, as pointed out 
by IGaskell & SnedderJ dl997l) . the ov erall sign i ficanc e 



is 



high. Early ex a mples in cluded [G askell (1988); 



transfer code ( Goosmann & Gaskel l2007p to model 
the effects on line profiles of scattering off an inflow- 
ing external medium. The two geometries considered 
are shown in Fig. 11. One is an infalling spherical dis- 
tribution of scatterers and the other an infalling cylin- 
drical distribution. In Fig. 12 we show a compari- 
son of observed profiles of two low- and high- ion- 
ization lines in PKS 0304-392 with various models 
with. We adopted an infall velocity of ~ 1000 km 
s" 1 based on velocity-dependent reverberation mapping, 
spect ropolarimetry, and the obs erved mean blueshift 
(see IGaskell & Goosmann! 12008 for details). It can be 
seen that both spherically and cylindrically symmetric 
models readily reproduce the blueshifting. 

An additional advantage of having significant scatter- 
ing in the BLR is that it sol ves the "smoothness prob- 
lem" for BLR line profiles ( Capriotti. Foltz. & B yard. 
Il98ll) . The intrinsic line broadening in an individ- 
ual BLR cloudlet is only of the order of the sound 
speed (~ 15 km s" 1 ), yet the velocity broadening of 
the BLR as a whole is hundreds of times greater. 
This requires the number of clouds to be ve r y hig h 



(Capriotti . Foltz. & BvardL 119811: lAtwood et all LL.982). 




The limit on the number of discrete clouds has now 
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&Gaskell| (|1989); |Crenshaw & BlackweU been pushed up to 10 8 dArav et allfl^bietrich et al 



Koratkar & Gaskell (1991a ffl; iKoristaetal 



bone & Kroli5 " 
More rece nt examples 



1996); 



can 



lUlrich & Horn. 
be found in 



Sere eev et all (11999b. iKollatschnvl (120031). iWelsh et all 
d2007l). boroshenko et al.1 d2008l). iBentz et al.1 (12008^ 



Den nev et al.1 (2009a) IBentz et all d2009cl) . Impor- 
tant independent evidence for inflow comes from 



high -resolution spectropolarimetry (e.g., ISmith et al 



2005). The systematic change in polarization as a 



function of velocity across the B aimer lines requires 
a net inflow of a scattering region somewhat exterior 
to th e B aimer lines. Pola rization reverberation map- 
ping dGaskell et al. . 2008a) can reveal the location of 
scattering regions. 

We have used the STOKES Monte Carlo radiative 



1 1999k . This constraint is relaxed if there is broadening 
by scattering. 

3.2. The overall velocity field of the BLR 

For a typical AGN, several independent lines of evi- 
dence (the blueshifting, velocity-resolved reverberation 
mapping, and spectropolarimetry) all point to the inflow 
velocity being of the order of ~ 1000 km s" 1 . As has 
been mentioned, velocity-resolved reverberation map- 
ping (see, for example, Fig. 9) implies that the domi- 
nant motion is not radial, but Keplerian or random. The 
observed widths of broad lines are indeed several times 



5 Available at http://www.stokes-program.info/ 




v (km s" ) 

Figure 12: Modelling the blueshifting of high-ionization lines. The 
profiles of OI A 1305 (narrow symmetric profile shown in red) and 
CIV ^1549 (thick blue line) for the quasar PKS 0304-392. The thin 
black line is the blueshifted profile produced by an infalling spheri- 
cal distribution of external scatters with an electron-scattering optical 
depth r es =0.5, and the dashed green line is the profile produced by 
the same distribution with r es = 1. The brown dots are the profile 
produced by a r es = 20 infalling external cylin drical distribution. The 
geome tries are shown in Fig. 11. Figure from Gaskell & Goosmann 
( 2008); PKS 0304-392 observations taken from Wilkes (1984). 



higher than the inflow velocity, and, of course, the pre- 
dominant motion for a flattened distribution must be Ke- 
plerian. 

As is clear from Fig. 6 and 7, when we observe 
the BLR (i.e., in type-1 objects) we are always seeing 
it close to face-on. The Keplerian component of ve- 
locity must be reduced by sin /, where i is the angle 
between the axis of rotational symmetry and the line 
of sight. The statist ics of line profiles in the SDSS 
dLa Mura et all l2009h suggest that for the vast major- 
ity of objects i is < 20 deg). 

As was realized by lOsterbrockl ( 1978k . the statis- 
tics of line widths imply that, in addition to Keplerian 
motion, there has to be a substantial additional com- 
ponent of velocity perpendicular to the orbital plane. 
Osterbrock appropriately called this "turbulence". The 
vertical component is also necessary for the recon- 
cilin g the structure of the BLR with i ts kinematics. 
As Mannucci. Salvati. & Strangal (1992) showed, for 
NGC 5548 the combined constraints of reverberation 
mapping and time-averaged line profile and favor the 
sort of "bird's nest" BLR distribution shown in Figs. 6 
and 7. 



locity components: 

VKepler > Vturb <L V» 



inflow 



(1) 



where the Keplerian velocity, v 'Kepler, of an emission line 
is a couple of times larger that the turbulent velocity, 
v tur b, which is in turn somewhat bigger than the inflow 
velocity, inflow. The ratios of BLR height to radius 
and of v Ke.n i e.r to V turb are similar to those deduced by 
Osterbroc 1<l <H978h . The only change to the Osterbrock 
model is recognizing that there is also a significant in- 
flow. 

4. Orientation effects 

It is well est ablished from radio properties (see 



Antonuccilll993l) that core-dominated AGNs are sim- 



ply lobe-dominated AGNs viewed from nea r the jet 
axis (i.e., near face-on). iGaskell et al.1 (120041) showed 
from a comparison of continuum shapes and line ratios 
that core-dominated and lobe-dominated AGNs have 
the same underlying optical-to-UV continuum shape 
and that the SED differences are just due to increased 
reddening in the lobe-dominated AGNs. We thus have 
every reason to expect the BLRs of core-dominated 
and lobe-dominated AGNs to be the same on aver- 
age. Lobe-dominated radio-loud AGNs should there- 
fore be an excellent laboratory for studyin g how orienta- 
tion a ffects the appearance of the BLR. iMiley & Milled 
1979) found that lobe-dominated AGNs preferen- 



tially had broader an d more irregular line profiles. 



Will s & Brown e (1986) discovered that the FWHM of 



HJ3 increases as we see AGNs more edge-on. This pro- 
vided strong support for a flattened BLR. 

AGNs with the peaks of their broad B aimer lines 
blueshifted or redshift ed from the systemic velocity 
have long been known (ILyndsl 1 19681) . It was proposed 
dGaskeli Il983l) that these peaks might represent sep- 
arate BLRs each associated with a member of a su- 
permassive black hole binary, but line profile variabil- 
ity observations on long and short timescales have de- 
livered two fatal blows to this hypothesis. Firstly, al- 
though for a while it looked like the expected binary 
orbital motion was being seen in long-term profile vari- 
ations in 3C 390.3 ( GaskellLll996l) . further observations 
showed that the radial velocity changes were completely 



inconsistent with a binary black hole dEracleous et al 



In summary, I believe that all the evidence points to 
the BLR having a nest-like appearance and having ve- 



1 19971) but were instead consistent with orbital mo- 
tion of concentrations of BLR gas orbiting in a disk. 
The second fatal blow was that velocity-resolved re- 
verberation mapping of 3C 390.3 strongly ruled out 
the binary BLR hypothesis because the redshifted and 
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Figure 13: The effect of broadening lines on the appearance of struc- 
ture in line profiles. The left frame shows a Lorentzian and two Gaus- 
sians chosen to approximate the appearance of Har or H/3 in 3C 390.3 
in 1981 or 1988. The right frame has the same line widths and peak 
intensities a s in the left frame, but ha lf the velocity displacements. 
Figure from Gaskell & Snedden (1997). 



blueshifted peaks varied simultaneously on a light- 
cross ing timescale (O'Brie net all 1 1 998L iDietrich et al 



1 19981) . This demonstrated conclusively that the double- 
peaked profiles arose from an inclined disk, as had 
be en widely proposed (see references and discussion 



Gaskell & Sneddenlfl999l) . Despite these double fa- 



tal blows to the idea that displaced broad-line peaks 
might be due to supermassive binaries, the topic of what 
signs there might be of sub-parsec supermassive bina- 
ries nonetheless remains one of considerable current in- 
terest (see review by Tamara Bogdanovic in these pro- 
ceedings). 

A subsequen t comprehensive survey of radio- 



galaxies by lEracleous & HalpernI ( 1994 (2003) revealed 
many disk-like B aimer line profiles. They found the 
FWHMs of the B aimer lines to be approximately dou- 
ble those of AGNs with single-peaked B aimer lines. As 
is shown in Fig. 14, a factor of two reduction in line 
width is sufficient to mak e displaced peaks disap pear. 



1999), IPopovic et ail d2004l) . and 



Gaskell & Snedden 

Bon et al] (120061) have argued that a disk-like emission 



line contribution is probably present in all BLRs but 
simply hard to recognize because, as illustrated in Fig. 
13, the classic double peaks become hard to see when 
the disk is near to face-on. 

It is straight forward to estimate the inclinations of 
the BLRs from broad disk-like line profiles. These 
can often be estim at ed to within a few degrees. 
Eracleous & HalpernI ( 1994 . 2003) get inclinations 



which predominantly have i > 25 deg. Their fits to the 
disk profiles also provide important confirmation that a 
significant turbulent velocity is needed and give the tur- 
bulent velocity for each object. Without the turbulent 
velocity component the peaks of the line profiles would 
be much too sharp. The turbulent velocities are fairly 



well determined from the line profile fits (to « ±250 
km s" 1 ). The average BLR turbulent velocity needed is 
1300 km s" 1 . This is roughly what would be expected 
from the height of the BLR/torus. The 1-cr scatter in the 
derived turbulent velocities is only ±400 km s" 1 , which 
is only slightly greater than the average formal uncer- 
ta inty i n the es timates. 

IBonI ( 20081) has estimated inclinations for single- 
peaked AGNs. For these we mostly see disks with incli- 
nations of i < 25 deg (see also Bon et al. in these pro- 
ceedings). The difference in sin i between the displaced- 
BLR-peak AGNs and single-peaked AGNs is thus about 
a factor of two. This agrees with the ratios of FWHMs 
for the two samples. 



5. The accuracy of AGN black hole mass determina- 
tions 

The component of velocity perpendicular to the equa- 
torial plane is vital for AGN black hole mass determi- 
nations! Without this strong turbulent velocity com- 
ponent, variations in sin/ would introduce substan- 
tial scatter into AGN black hole mass estimates, espe- 
cially since type-1 AGNs are observed close to face- 
on. There is recent evidence that there is remarkably 
little scatter in AGN mass estimates. Firstly , it has 
become apparent (Bo chkarev & Gaskell 20091) that the 
two main methods of estimating black hole masses from 
the BLR agree surprisingl y we l l. Th e Dibai single- 



1977 ) and reverbera- 



ly we i 

epoch- spectrum method (iDibail 
tion map ping methods agree to within the expected er- 
rors. Gaskell (2009b) has shown furthermore that a 



simple refinement of the method produces even better 
agreement. The agreements mean that such methods are 
estimating the effective rad ii of the BLR correctly. As 
Bochkarev & Gaskelll (2009) discuss, the success of the 
Dibai method means that the inner regions of AGNs are 
very similar. In particular: 

1. The spectral energy distribution (SED) from the 
optical to the far UV must be very similar in all 
type-1 AGNs because the optical region where the 
flux is measured is far removed in energy from 
the far UV which is photoionizing the gas. Al- 
though A GN SEDs look different iGaskell et al j 
d2004l) and lGaskell&Benkerl d2007l) have already 
argued that the apparent variation is not real but is 
primarily caused by reddening. 

2. There is a simple scaling relationship between 
the luminosity and the effective radius. This is 
supported by reverberation mapping estimates of 
the effective radii of BLRs dKoratkar & Gaskell. 



199ld:lKaspi et all l2000L [20051: bentz e7aill2006l 



2009a) 



Both the Dibai and reverberation-mapping methods 
of estimating black hole masses depend on observed 
BLR line widths, so geometric differences and orien- 
tation effects will affect both methods. An important 
external check on the accuracy of AGN black hole mass 
estimates is provided by the tightness of the relationship 
between black hole mass, M., and luminosity, Lhost, 
of th e bulge of the host galaxy. iGaskell & Kormendy 
(2009) have recently shown that estimating M. by the 
Dibai method and Lh os t from the fraction of starlight in 
SDSS spe ctra gives a scatter of ±0.23 dex in log M. (see 
Fig. 14). iBentz et al.1 (I2009bl) have estimated Lh 0S t com- 
pletely independently for a different set of AGNs using 
HST photometry and published reverberation mapping 
mass estimates. They get a scatter in logM. of +0.33 
dex. Both of these scatters in the AGN M. - \ogLh os t 
relationships are sm aller than the +0.38 dex scatter 
iGultekin et al.1 d2009h and others find when M. is deter- 
mined by stellar dynamical methods, but they are still 
greater than the +0.17 dex scatter in the M. -cr* rela- 
tionsh ip for pure bulge (i.e., barless) galaxies dGrahami 
2008). T he Dib ai method and the method proposed by 



Gaskell ( 2009b) seems to give particularly tight M. - 
<x* and M. - Lhuig e relationships f or the most massive 



elliptical galaxies ( GaskellL 2009ak . This is probably be- 
cause they have the least intrinsic scatter in the M. - <x* 
relationship. These comparisons with predictions from 
host galaxy properties imply that black hole mass de- 
terminations from the BLR are surprisingly accurate - 
as accurate as the best stellar-dynamical estimates. This 
accuracy of black hole mass estimates made using the 
BLR provides strong support for all type-1 AGNs be- 
ing very similar as far as the structure and kinematics 
of the BLR goes, and for orientation effects being mini- 
mal. The accuracy of AGN black hole mass estimates is 
thus consistent with there being a substantial turbulent 
BLR velocity component and type-1 AGNs being seen 
close to pole-on. 

6. What drives BLR motions? 

The circular component of motion of BLR clouds 
is a simple consequence of gravity, but the turbu- 
lent vertical component of motion has to be main- 
tained against dissipation losses, and an outward trans- 
fer of angular momentum is necessary to get in- 
flow. We now know that the viscosity needed to drive 
the outward flow of angular momentum in accretion 
discs, and hence the inward flow of matter, is the 
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Figure 14: T he M. - L host relationship fo r 100 AGNs with 0.13 < 
z < 0.34. See lGaskell & Kormendvl d2009h for details. M. has been 
estimated by the Dibai method. The diagonal line is the OLS-bisector 
fitM.ocL?- 8 ! 



magneto-rotation instability (MRI) (Balbus & Hawley, 
Il99ll) . Over the last decade increases in comput- 
ing power and the development of more sophisti- 
cated programs by several groups have allowed increas- 
ingly detailed magneto-hydrod ynamic (MHD) simula- 

2001; 
20051 



tions of accretion flows (e.g., Hawlev & Krolik 

1 1 I — — I I ^ — — 

Prosa 2003; Anninos, Fragile, & Salmonson 

lOhsuga et al.l 120091: IShafee et al.l 120081 and references 
therein). In these models attention has been focused 
a lot on the low-density outflows. Because emissivity 
goes as the square of density, the emission is dominated 
by the high-density. To my mind what is impressive 
about every single one of these models is that, despite 
the different modeling approaches, the velocity fields of 
the high-density material all match the velocity field in- 
ferred for the BLR. I.e., the dominant motion is Keple- 
rian, but there is substantial turbulence, and a significant 
inflow. This is very clear when one watches movies dif- 
ferent groups make of their simulations. I believe that 
these simulations give us a physical basis for what we 
have deduced from BLR observations: the BLR is the 
material accreting onto the black holes. It was indeed 
noted a long time ago that if the BLR is inflowing it can 
provide the necessary mass flux for powering the AGN 
(Pad ovani&R afanelli. 198§). 



7. Conclusions and unsolved problems 

I think we now have a fairly good emerging picture of 
what the BLR is like and what role it plays in the life of 
an AGN. Interstellar material approaching the nucleus 



settles into a flattened distribution, the thick torus. Ma- 
terial loses angular momentum because of MRI turbu- 
lence and gradually spirals inwards. When the mate- 
rial of the torus gets within the dust sublimation radius, 
the dust evaporates and we have the BLR. The turbulent 
BLR continues to spiral inwards towards the black hole 
where it is eventually accreted. The degree of ionization 
increases as the gas gets closer in. The optically-thick 
material, which will tend to be concentrated towards 
the mid-plane, produces continuum emission; the more 
optically-thin material produces the BLR. Not all of the 
BLR is accreted. Some of it is driven off the surface of 
the BLR/torus in a high- velocity, low-density wind, as 
is found in all the MHD simulations and as is observed. 

Although I believe we are getting a clear overall pic- 
ture of the BLR, there is still plenty to work on both ob- 
servationally and theoretically! For a start, the picture 
discussed above needs to be thoroughly tested to ver- 
ify that it works for all objects and not just a few well- 
observed objects such as NGC 5548. More work needs 
to be done to see whether a disk- wind model (the lead- 
ing rival to the model presented here) could also explain 
everything. We know that there is outflowing gas as 
well as gas accreting onto the black hole. The question 
is: how much of this is also contributing to the broad- 
line profiles? ( especially to the high-ionization lines). 
Illic et al.1 d2008l) . for example, have shown that outflows 
can match some observed BLR profiles, so determin- 
ing the relative contribution of an outflow to broad- line 
profiles from line-profile fitting alone is difficult. I think 
that reverbera tion mapping and spectropolarimetry (see, 
for example. lAxon et al.ll200 8) are going to provide the 
best answers. Kollatschnv (2003) found marginal evi- 



dence f or some outflow in Mrk 110, and iDenney et al 



(2009b) have found a clearer signature of an apparent 
outflow component in velocity- resolved reverberation 
mapping of NGC 3227. Since the bennev eTaD d2009bl) 
results are from a single short observing campaign, I do 
not think that NGC 3227 presents a major problem yet 
for the general picture presented here@ If follow-up ob- 
serving campaigns confirm the signature of outflow in 
NGC 3227 then this would be a significan t challenge to 



the model favored here. Nevertheless, the Pe nney et al 



2009b) result does caution us that AGNs might not all 



6 The uncertainties in the red-wind/blue-wing lags can be larger 
than thought. NGC 5548 provides a good illustration of this. The red- 
wing/b lue-wing lag varies from year to year by more than the formal 
errors dWelsh et all l2007h . but the NGC 5548 BLR is probably not 
changing direction at the end of every observing season! A strong 
reason for believ ing that the N GC 3227 kinematics are not unusual is 
that, as Dennev et al. ( 2009b) point out, the mass estimate lies on the 
M. - <t* relationship. 



be identical in the relative dominance of inflow and out- 
flow. 

Even if we are right about the basic structure of the 
BLR and torus of AGNs, there are still a lot of inter- 
esting and potentially important details in need of fur- 
ther investigation. Although in this review I have been 
emphasizing the similarities among AGNs and what 
they imply, there are some significan t differences i n the 



BLRs too (see, for example, Marziani et al. 19961) . If 



our basic framework of how an AGN works is correct, 
then the differences need to be explicable within the 
framework too. Space here only permits a brief men- 
tion of some of these problems, but fortunately many 
of them are reviewed and discussed elsewhere in these 
proceedings (see, for example, the reviews by Mike Er- 
acleous and Jack Sulentic). 

It has been known for over three decades now that 
object-to-object differences are correlated with each 
other, and one of main drivers of the correlated d iffer- 
ences is the Eddington ratio (see Sulentic et al.l2000l and 
these proceedings). Since we now have reliable AGN 
black holes masses, we also have reliable Eddington ra- 
tios, so there is a lot that can be done in investigating the 
dependence of BLR properties on accretion rate. I think 
there is a lot that needs explaining here. 

The biggest object-to-object d ifference in optica l 
spectra is optical Fell emission dOsterbrockl 1977 ). 
Understanding how the very strong optical Fe II emis- 
sion seen in AGNs i s produced has been a long- 
stand i ng problem (se e iBaldwin et al] 120041: Ijolv et al 
20081; iHu et al j 120081: iKuehn et al j 120081: IVerner et al 
20091: iDong et al.ll2009l for recent discussions). In the 



BLR model discussed here, optical Fe II emission arises 
in the outer part of the BLR just inside the torus (and 
quite likely overlapping with it), but this does not read- 
ily explain why optical Fe II is so much stronger in some 
objects than others. 

Another mystery of the correlated object-to-object 
differences is the strength of the narrow-lint region 
(NLR) emission. This is the other strong object- 
to-object difference and it is mys teriously st rongly 
anti-correlated with Fell emiss i on dOsterbrockl 1 1977 1 : 
Steiner, 1981 : boroson & OkeL Il984l: iGaskelll. 1 19871: 



Boroson & GreenL 119921) . A complete model of AGNs 



needs to explain why the NLR and BLR know about 
each other. 

Although I have argued that the basic properties of 
AGNs with broad disk-like B aimer line profiles are con- 
sistent with the picture presented here, these objects, 
and especially the variability of their profiles, present 
some special challenges, as is discussed in Mike Era- 
cleous's review. There is also a lot more to be learned 



10 



with orientation effects. 

In summary, I think that although our overall pic- 
ture of the BLR and the role it plays in the AGN phe- 
nomenon is becoming clearer, many mysteries remain, 
there is still a lot to learn, and there are probably sur- 
prises in store. 

I am grateful to Luka Popovic, Milan Dimitrijevic, and 
the other members of the scientific organizing com- 
mittee of the 7th Serbian Conference on Spectral Line 
Shapes for inviting me to speak on this topic. I would 
like to express my appreciation to Dragana Ilic and 
all the members of the local organizing committee for 
providing a very pleasant and stimulating experience 
throughout the conference, both culturally and scientif- 
ically. I also have to thank my collaborators and for- 
mer graduate students for all their contributions and dis- 
cussions over the years, and the anonymous referee for 
useful comments. This research has been supported in 
part by US National Science Foundation grant AST 08- 
03883. 
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